Liquid meniscus friction on a wet plate: Bubbles, lamellae and foams 
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Many microfluidics devices, coating processes or diphasic flows involve the motion of a liquid 
meniscus on a wet wall. This motion induces a specific viscous force, that exhibits a non-linear 
dependency in the meniscus velocity. We propose a review of the theoretical and experimental work 
made on this viscous force, for simple interfacial properties. The interface is indeed assumed either 
perfectly compressible (mobile interface) or perfectly incompressible (rigid interface). We show that, 
in the second case, the viscous force exerted by the wall on the meniscus is a combination of two 
power laws, scaling like Ca 1 ^ 3 and Ca 2 ' 3 , with Ca the capillary number. We provide a prediction 
for the stress exerted on a foam sliding on a wet solid and compare it with experimental data, for 
the incompressible case. 

Many industrial or natural situations involve the motion of liquid menisci on a wet solid. For enhanced oil recovery 
[IHIj], or soil remediation liquid foams are pushed into porous media to displace oil or pollutants; surfactant 
lamellae are produced in the lung and move in the bronchial tubes 0, Q ; innovative set-ups are developed to control 
bubbles motion in micro-channels for lab-on-a-chip applications Q ; pulling a solid out of a liquid bath is a common 
way used in industry for coating [8j. 

As depicted in Fig. [TJ the same local flow is observed close to the wall in all these examples. While the meniscus 
moves on the wall, it deposits a thin wetting film. Its thickness, of the order of few tens of microns, depends on 
the meniscus velocity through the capillary number Ca = r\U /j that compares the viscous stress and the capillary 
pressure, through the values of the solution viscosity ry, the surface tension 7 and the meniscus velocity U. This 
property, that makes possible a finely tunable coating, has been widely investigated [8]. In contrast, the force required 
to pull a meniscus at a given velocity on a wall is much less understood, and plays a crucial role for the motion of 
confined bubbles or confined foams. This force is governed by the velocity gradients in the liquid phase. As the liquid 
phase is confined into the small gaps between the bubbles and the wall, these gradients, and consequently the viscous 
force, are strongly enhanced and depend on the dynamical shapes of the meniscus and of the wetting film. As a direct 
consequence, pushing a foam in a tube requires a pressure drop that is much higher than the one needed to push the 
surfactant solution in the same tube at the same mean velocity @,[l0|. 

The theoretical prediction of the shape of the deposited film has been proposed by Landau and Levich (TlJ and 
by Derjaguin [12] (LLD) in the limit of small capillary numbers. Brctherton [l3[ and Park and Homsy [l4| adapted 
this calculation to the case of a confined bubble and predicted the pressure drop associated with the bubble motion. 
Experimental results mainly report on the film thickness deposited by long bubbles in cylindrical tubes [H, Il5l - [20j or 
on solids pulled out of a bath (Q and references therein, 21]), for a typical range of capillary number [10~ 6 ; 10 -3 ]. 
The viscous force acting on the meniscus, or the pressure drop, have been measured for different systems including 
bubbles separated by liquid slugs [22|, [23[ or lamella [ll [23 - |27| moving in tubes, and foams in 2D [28| or 3D geometries 
[29l |30| . Numerical simulations have extended the pre dictions for the meniscus induced pressure drop to higher 
capillary numbers or to more complex geometries [3ll - l39l ] . 
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FIG. 2: Schematic view of a liquid meniscus on a solid that moves in the :r-direction at velocity U. The thickness profile is 
h(x), the pressure is p(x) and the x-component of the velocity in the meniscus frame is v(x,y). The radius of curvature of the 
meniscus tends to r m for x — > — oo (meniscus side) and the film thickness tends to /itx> for x — » oo (wetting film side). The 
plate velocity U > corresponds to a front meniscus (f) and U < to a rear meniscus (r). The problem is invariant in the z 
direction. 



FIG. 1: Different diphasic flows involving the motion of a meniscus on a wet wall. From left to right: a single bubble in 
a tube (Bretherton's problem [13]) or in a Hele-Shaw cell; a lamella in a tube; a plate or a fiber [8J pulled out of a liquid 
bath (Landau-Levich-Derjaguin (LLD) problem [TH |4C(|): a droplet spreading on a wet plate (Tanner's problem [41]). (7 is the 
velocity of the solid in the frame of the meniscus. The front and rear menisci are denoted respectively by (f) and (r). Each 
meniscus is connected to a wetting film of length I. 

Lamella and foams are stabilized by surfactants that induce a specific stress at the interface. If the area A of an 
interface element increases fast enough, the surfactant concentration decreases and the surface tension increases, thus 
acting against the extension. The dilatational modulus of the interface, defined as E = (l/A)dj/dA, is a measure 
of this phenomenon. Moreover the intrinsic viscosity rj s of the surfactant monolayer induces an additional interfacial 
stress, that also acts against the extension. In the limit of high dilatational modulus or high interfacial viscosity, the 
interface behaves as an incompressible surface. 

In this paper, we discuss the motion of a meniscus on a wet wall, in the ideal cases of incompressible interfaces or 
interfaces with uniform tangential stress. The results obtained in these limiting cases constitute important references 
for models with more realistic interfacial properties. SectionUpresents classical calculations for the dynamical meniscus 
shape, from which the viscous forces are deduced in Section HH as well as the validity range of the ideal models. An 
original derivation of the pressure drop for a single bubble is proposed in Section Hill and compared with Brctherthon's 
calculation. Finally, adapting Dcnkov's approach [29(, these results are used in Section HVl to derive predictions for 
the foam/ wall tangential stress as a function of the bubble size and liquid fraction, for stress- free or incompressible 
interfaces. For incompressible interfaces, we show that the stress is not a power law of the capillary number as usually 
assumed, but rather a sum of two contributions scaling as Ca 1 / 3 and Ca 2//3 , in good agreement with experimental 
results. 



I. DYNAMICAL MENISCUS SHAPE 
A. Resolution of the lubrication problem 

We consider a liquid meniscus moving at constant velocity in the x direction on a solid plate, in steady state. We 
assume an invariance in the z direction (see the convention of orientation on Fig. [2| and we neglect inertial and 
gravitational effects. In the meniscus frame, a positive plate velocity U > corresponds to a plate pulled out of a 
bath as in the Landau Levich problem, to the front end of the bubble as in the Bretherton problem or to the rear side 
of lamella (see Fig. [I}. Keeping the Bretherton notation, this case is called the front meniscus, denoted by / in the 
following. Conversely, U < corresponds to a plate plunging into a bath, to the rear end of the Bretherton's bubble, 
to a drop spreading on a wet plate or to the front side of a lamella, and is called the rear meniscus, denoted by r. 

This problem is governed by the Stokes equation, in the lubrication regime. We denote the pressure in the liquid 
phase p(x), the velocity in the x direction v(x, y) and the film thickness h(x). The fluid viscosity is rj and the surface 
tension 7. The equation of motion can then be written as: 

d 2 v I dp 7 d 3 h 
dy 2 rj dx 77 dx 3 
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Using v(x, 0) = U this leads to 



The continuity of tangential stress er mt provides the boundary condition at the interface. 

da int dv 

The exact expression for a mt depends on the physico-chemical properties of the solution. This interfacial stress 
is denoted by 7 in eq. [1] which implicitly assumes that, at dominant order in Ca, 7 = j eq , the equilibrium surface 
tension. The explicit notation 7 eg and a mt will be used to discuss the small difference between both terms. Elsewhere 
we will keep the generic notation 7. Without surfactant <r mt = j eq is exactly verified, and the relation ^ imposes 
= at the interface. This case is the stress free case denoted (sf) in the following. In the opposite limit of 

incompressible interface, the boundary condition is div 2D u = 0, with div 2D the divergence in the plane tangent to 
the interface and u the interfacial velocity. Consequently, small variations of the interfacial stress occur along the 
interface, discussed in section III B 21 In our geometry, because the flow is quasi-parallel and invariant in the z direction, 
&\v 2D u simplifies into dv(x,h(x))/dx. The incompressibility of the interface thus implies that v(x,h(x)) is constant. 
This uniform velocity v(h) is imposed by an external condition, in the static meniscus or in the wetting film and 
may a priori take any value, depending on the whole experimental situation. In this paper, two particular values are 
considered. The case v(h) = v(0) = U will be called the rolling case (rol): in the region of interest, close to the wall 
and to the meniscus, the interface makes a rotation around the bubble in the bubble frame and has no relative motion 
with respect to the wall. The case v(h) = is the sliding case (sli): the interface is locally immobile in the bubble 
frame and slides on the moving wall. 

In order to take into account these three different boundary conditions (sf, sli and rol) within the same formalism, 
we use a general velocity expression that depends on the two parameters A and /z, as summarized below: 



v=^(\hy-y 2 /2) + U(l-w/h). (4) 

Finally, we impose the condition that the wetting film is flat far from the meniscus, with a thickness hoa- The flux 
(per unit length in the z direction) in the flat region is Q = f7/ioo(l — fi/2) and flux conservation leads to the central 
relation 

8 3 h 

± 3/3Ca(^oc - h) = -^h 3 (5) 

with /3 = (2 — /i)/(3A — 1) and Ca = rj\U\/j. By convention, we use the upper sign (here +) for the front meniscus 
and the lower sign (here — ) for the rear meniscus. 
For the different boundary conditions we obtain: 

• Stress free case (sf): A = l, (J, = Q, (8 = 1 

• Sliding case (sli): A = 1/2, p = 1, /3 = 2 

• Rolling case (rol): A = 1/2, fi = 0, /? = 4 

We now define the dimensionless quantities h = hooH, x = h oc X/(3(3Ca) 1 ^ 3 , and rescale equation (JSJ into 

H"'H 3 = ±(1 - H) (6) 

The boundary conditions are H(oo) = 1, H'(oo) = and H"(oo) = 0. The problem is finally closed with the matching 
condition h"(—oc) — l/r m , with r m the static meniscus radius of curvature. This condition is equivalent to 

hoo =r m (ZpCa) 2 / 3 H" (-00). (7) 

A numerical resolution requires the determination of boundary conditions at finite distances, and these are obtained 
from the linearized problem. At large X, H is close to 1 and the problem can be linearized using a small parameter 
e(X) = H(X) — 1. The governing equation © becomes e'" = =Fe, which has solutions of the form e = e KX , with 
K 3 = =Fl. For the front meniscus K = — 1 is the only root with a negative real part and the solution H(X) = l+eoe~ x 
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can be used to determine the values H(0), H'(0) and H"(0) needed to solve eq. [6] numerically. Changing the arbitrary 
value of £o only translates the solution along the X direction. The unique solution (X) is shown on Fig. [3] (inset). 
The numerical value of H'l(—oo) is H'J t = 0.6430, and, from eq. [3 we recover the classical results [H} 

h^ = 1.34 r m Ca 2 / 3 ; h%=2 2 / 3 h°J > ; Kg = 2 4 / 3 ^ . (8) 

In the case o f a plate pulled out of a liquid, the result is usually expressed as a function of the capillary length 
= v/7 '/{P9)i with p the solution density and g the gravity. Due to the hydrostatic pressure, the meniscus curvature 
changes with the height and can be shown to be r m — l c /^/2 at the top of the meniscus. The film thickness is thus, 
in the stress free case, = 1.34Z c Ca 2 / 3 /\/2 = 0.94 Z c Ca 2 / 3 . 




FIG. 3: Family of rear meniscus profiles, obtained by solving eq. ([6]). The boundary conditions at X — are derived from 
eq. [9] using eo = 10 -6 and different values of <6. The dashed line (red online) is the r* solution, having the same asymptotic 
curvature as the front meniscus. Inset: solution of eq. ^ for the front meniscus. 

For the rear meniscus, the linearized equation e'" — e admits two roots with a negative real part K = — l/2±iv3/2 
. The general solution of the linearized equation is thus 

e = e a e- x/2 cos(\/3X/2 + $) . (9) 

Different values of Sq just translate the solution, but different values of $ lead to different solutions, with H" converging 
to a large range of possible values (positive or negative) when X —> — oo 43]. The family of rear meniscus profiles are 
represented on Fig. |3] This states that, for a given meniscus radius of curvature, we can find a steady solution for 
any film thickness initially present on the wall. In the following discussion, the profile with H"(—oo) = 0.643, as for 
the front meniscus, will be denoted by the subscript r*. It plays an important role, as it is the shape obtained when 
the rear meniscus moves on the wetting film deposited on the wall by a front meniscus of same radius, moving at the 
same velocity. 



B. Asymptotic parabolic shapes 

The asymptotic expression for the film thickness for X — > — oo is a parabolic profile H gu t(X) = H"{—oo){X — 
X p ) 2 /2 + Hp, with X p the position of the minimum of the parabola and H p — H out (X p ) [44j] . Both X p and H p need 
to be determined numerically. The best numerical fit, plotted on Fig. |H have been obtained with = 2.88 for the 
front meniscus with X in the range [—400;— 150]. A slightly smaller value, H^ Br = 2.79 is reported in [13j and is 
found when the fit is performed for X in the range [—90; —40]. 

For the rear meniscus, we found = —0.82. Its asymptotic profile thus intersects the wall at the contact point X c 
shown in Fig. 2J The apparent contact angle 6 r al also shown in Fig. 21 has been computed for any initial wetting film 

thickness in (4{|. It can be deduced from the parabola equation: H' r out (X c ) = 2 H'^ out and from the relation 
{ipCa) 1 ^ H' r out {X c ) — tan 9 r a ~ r a . For the three boundary conditions considered, we get, for the r* solution: 



9 s J' r ' =lACa 1/3 ; ef' r ' =1.7Ca 1 / 3 ; Q r °^ r ' = 2.2 Ca 1/3 (10) 
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These relations are ver y si milar to Tanner's law, which describes the spreading of a droplet on a solid wall in the 
perfect wetting case |4ll l46l - l48j . However, for a spreading droplet, the profile far from the wall is not parabolic, 
but rather close to a straight wedge. In that case, the apparent contact angle involves a logarithmic term. Indeed, 
imposing the boundary conditions = when x — > oo and h" (x) — > when x — > — oo in cq. ([5]) leads to [ill . l46j 

9 d a rop (x) = [^Hx/x^Y^Ca 1 ^ . (11) 

The micro-scale cut off is obtained by matching the film profile to a precursor film initially present on the wall and 
governed by Van der Waals forces : x m i n = 2.5A/(QjT^/)3^ 1 ^ 6 Ca^ 2 / 3 , with A the Hamaker constant [46j . 




FIG. 4: Asymptotic parabolic shapes of the meniscus. The solid lines are the numerical solutions of eq. [6]for the rear meniscus 
r* (left) and for the front meniscus (right). The dashed lines are the best parabolic fits. The horizontal dotted lines show the 
height of the parabola minimum Hp — —0.82 and H^ = 2.88, for the rear and front cases, at the position X p . The apparent 
contact point of the rear meniscus is X c and the apparent contact angle at this point is Q r a . 



C. Comparison with some experimental data 

1. Wetting film thickness in the LLD geometry 

The wetting film thickness as a function of the plate velocity has been extensively studied for the front meniscus in 
the LLD geometry, with fibers or plates pulled out of pure liquids or surfactant solutions. The stress free result is 
obtained for pure liquids. For surfactant solutions, the whole range of thicknesses between and h™ 1 are observed, 
depending on the surfactant concentration, the imposed velocity and, for the fiber case, the fiber radius [1, Hlj]. In 
LLD geometries, surfactant solutions with high surface modulus lead to a wetting film thickness in good agreement 
with the rolling case. While the solid is pulled out of the liquid reservoir, the total interfacial area increases: the 
new area is produced at the free interface of the liquid reservoir, in agreement with the rolling assumption. This is 
confirmed by recent velocity field measurements in the reservoir [50j . 

The sliding case implies that the whole wetting film is sheared. It is thus incompatible with a very long wetting 
film (see section III B 31 for a quantitative discussion) . It is therefore not observed in the LLD geometry, but can be 
obtained for bubbles, as discussed below. 



2. Film profile for bubbles with incompressible interfaces 



The film profiles shown in Fig. [5] were measured by Denkov et al. using interferometry, for a solution of very high 
dilatational surface modulus [51| . We compare these experimental results with the rear and front solutions of eq. ([6]) . 
As the meniscus curvature r m is the same on both sides we use the r* solution for the rear meniscus. The coordinates 
x = hooX/ (SflCa) 1 / 3, + xo and h = h^H are deduced from the numerical solution H(X) using the experimental 
capillary number Ca, the film thickness (measured in the middle of the film) and a translation parameter xq that 
is independently fitted for the rear and front parts of the profile. The three possible values of j3 have been tested. For 
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the two slower cases, the rear and front solutions nicely overlap on the central region and a very good fit is obtained 
using /3 — 2 (sliding case). The profiles obtained for (3 = 1 (stress free) and /3 = 4 (rolling) are also plotted on Fig. 
EJinset) for one meniscus and show a significatively worse agreement with the experimental profile. It thus confirms 
that the foaming solution used, containing insoluble surfactants, leads to a sliding regime in this geometry, as already 
assumed in [5l| . 

At larger velocities, the contact area between the bubble and the wall, i.e. the area of the central flat film, decreases 
and eventually disappears. The two solutions H r (X) and (X) are both determined by assuming a flat film far 
from the meniscus and are thus not valid anymore. In that regime, the whole hydrodynamical problem, involving the 
two meniscus at the same tim* 3 sbnnlH he anlwd 





FIG. 5: Wetting film profiles for a 2D bubble in contact with a wall moving to the right, for incompressible interfaces. +, x 
and •: experimental data from Fig 12a in [Hlj], for Ca = 1.15 10 -4 , 2.3 10 -4 and 4.6 10 -4 , respectively. Full lines: solutions of 
eq. (|6} (shown in Fig. [3| rescaled using /3 = 2 (sliding case). Insert: Zoom on the rear meniscus of the first profile. Dashed, 
full and dash-dotted lines: solutions of eq. (|6j rescaled using j3 = 1 (stress free case, red online), /3 = 2 (sliding case) and f3 = 4 
(rolling case, blue online) respectively. 



3. Flux around a bubble for pure water 



An indirect measurement of the wetting film thickness has also been obtained by Schwartz, Princen and Kiss 18] 
for the case of a bubble in pure water. They pushed a single bubble along a tube of radius r t and measured the 
solution flux in the bubble frame Qtot- They report a fractional speed difference W defined by W — Qtot/ '(nr^U). 
From eq. gj we predict Q to t = 2-Kr t Uhoo{l - /i/2) and W = 2h 00 (l - (i/2)/r t . Finally we get W sf = 2.68Ca 2 / 3 , 
W sh = 2~ 1 / 3 W s f and W roL = 2 4 / 3 W s f , that we compare with the experimental data in Fig. |H1 The flux measured 
for short bubbles is in nice agreement with the stress free prediction, as expected for pure water. However, for long 
bubbles (more than 1.5 cm for a 1mm inner diameter tube) the flux increases progressively, which is in contradiction 
with a transition toward a sliding case, for which a smaller flux is expected. In that case, the experimental observation 
is instead compatible with a transition toward a rolling observed in LLD geometry [521 ] . 

It is not clear why the regime is modified with the bubble length in this case and, more generally, what determines 
the sliding or rolling regime in the limit of incompressible interfaces. Note that for a bubble in a cylindrical tube, the 
flow is necessary axisymmetric. The rolling case is thus not compatible with the incompressibility of the interface at 
the global scale, as some interface needs to be created at the front and suppressed at the rear of the bubble. 



bubble length (cm) 



FIG. 6: Fractional speed difference W of the liquid phase around an isolated bubble in a tube. • : experimental data from Fig. 
4 with Ca — 3 10~ 5 and rt = 0.5mm. The solid lines are the predictions, from top to bottom, for the rolling case, the 
stress-free case and the sliding case. 



II. VISCOUS FORCES IN A MENISCUS 



A. Theoretical predictions 

The viscous force exerted by the solid wall on the liquid (per unit length of meniscus in the z direction) can be 
deduced from the velocity expression eq. [4] 

f f/r = -V ^(x,0)dx=-ri \1 hdx + vUv -rdx. (12) 

J meniscus J — oo V ®X J —oo " 

The first term converges at +oo and — oo and corresponds to a viscous stress localized in a small domain close to the 
apparent contact point between the meniscus and the wetting film, as shown in dimensionless form on Fig. [7J This 
domain is usually called the dynamical meniscus, and its extension Id scales with h^/ (S^Ca) 1 ^ , i.e. as rmCa 1 ^, 
with a prefactor discussed in section Hi B 31 

The second term diverges at +oo and only appears for the sliding case (fi = 1) for which the whole wetting film is 
sheared. It thus imposes a cut-off length I corresponding to the actual length of the wetting film. Expressed in term 
of the universal function H the force (per unit length) is given by: 



ff/ r = 7 (3/3Ca) 2/3 



\F f ' r + J±-G f/r (L) 
op 



(13) 



with Ff/ r = =f 1^(1 - H)/H 2 , G f ' r {L) = ± l/H and L = (£/r m )(3/3Ca)" 1 / 3 / 'H'^ out . As L depends on Ca, 
the last term introduces some additional Ca dependence into the viscous force. 



In the front meniscus, H > 1 everywhere, implying that the viscous force exerted by the wall is positive, as expected 
for U > 0. The numerical integration leads to = 1.85. In comparison, the rear meniscus is oscillating around 
H = 1 and, as underlined in [13], the sign of the force depends on the value of H"(—oo). We discuss this point in 
section ITlI A II The value of F r is plotted as a function of H"(—oo) in Fig. [51 For the r* profile, we find F r * = —0.53 
and thus, as for the front meniscus, a viscous force exerted by the wall in the direction of its own motion. The 
vanishing force corresponds to H"(—oo) — 1.1, which implies hoo = 1.1 r m (3/3Ca) 2 / 3 . 

The wetting film contribution G^ r (L) tends to an affine function of L at large L, because H(X) becomes constant. 
The coefficients of the affine law depend on the precise definition of L. In this case, we choose L as the distance 
between X p , the minimum of the asymptotic parabolic profile (see Fig. 2]), and the end of the wetting film. With this 
definition, we get, for long wetting films: 

G f — L + 1.27 and G r * = -L - 4.0 (14) 

The agreement with the exact numerical solution is within a few percent for L > 2. 

Finally, the viscous force per unit length for each investigated case is, in physical units: 
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FIG. 7: Normal and tangential force distribution in the front and r* meniscus, in dimensionless units, (a) Upper graph: Front 
meniscus profile. Lower graph: solid line: pressure P — 0.643 — H"(X); dashed line: meniscus contribution to the viscous 
stress (H — 1)/H 2 ; dot-dashed line: velocity divergence —H'/H 2 . (b) Same functions for the rear meniscus r*. Full line: 
P = 0.643 - H"(X); dashed line: (1 



12 3 4 5 

H"(-inf) 



FIG. 8: F r — J^ a (x> (l — H)/H 2 dX as a function of the asymptotic curvature H"(—oo). F r is the dynamical meniscus contribution 
to the viscous force per unit length exerted on the rear meniscus. It is related to the corresponding physical quantity by the 
eq. 1131 The circle represents the r* solution. 



ff f = 3.84 7 Ca 2 / 3 /;; = -1.1 7 Ca 2 / 3 (15) 

ff u = 3.75 7 Ca 2/3 + 0.47 7 Ca 1/3 — f r M = -3.07 7 Ca 2 / 3 - 0.47 7 Ca 1 / 3 — (16) 

f f rol = 4.84 7 Ca 2 / 3 ff ol = -1.38 7 Ca 2 / 3 (17) 

Given the orientation conventions used in the paper, for a wall moving in the positive direction, the force per 
unit length exerted on a bubble or on a lamella by the wall is given by / = ff — f r . We define the parameters 
£ = (4.94; 6.82; 6.22) and £ = (0; 0.94; 0) respectively for the stress free, sliding and rolling cases, so that 

f = ^Ca 2 / 3 + ^Ca 1 / 3 — . (18) 

The viscous force exerted on a meniscus thus exhibits a power law in relation to the velocity, with an exponent 
of 2/3 for the stress free case and for the rolling case. However, for the sliding case, the force is a combination of 
two power laws of exponents 1/3 and 2/3. For a typical capillary number Ca = 10~ 3 , both terms are of the same 
order of magnitude if the meniscus size r m is 10 times smaller than the bubble size I. This latter result provides 
some theoretical clarification to the seminal work by Denkov et al. which first shed light on the specific behavior of 
foams with incompressible interfaces. The approximate resolution of eq. [5] for the sliding case developed in [29J relies 
on an approximated shape for the wetting film and on the assumption that its thickness is governed by a minimum 
dissipation principle. It leads to the prediction that the contribution of the sheared wetting film to the friction force 
scales as 7 (^/r m ) 1 / 2 Ca 1 / 2 . A 1/2 power law is found for deformable objects with elastic properties but, as shown 
in this paper, it does not apply to the friction between bubbles and a solid, even in the incompressible interface limit. 

The two power laws appearing in eq. (1181) can be derived from simple scaling arguments. The term localized in the 
dynamical meniscus is built from a velocity gradient U/h^, integrated along the length id ~ h oc Ca^ 1 /' i . This leads 
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to a force per unit length varying as rjUCar 1 / 3, — ^Ca 2 ! 3 . The term obtained from the shearing of the whole wetting 
film of length £ scales as rjiU/h,^ ~ ^(l/r^Ca 1 ^ 3 . 

Finally, the viscous forces obtained in the rolling and stress free cases are simply related to each other by the 
relation f^jj — 2 1 / 3 f^j r . This is not the same rescaling as for the film thickness h^, found to be /i™' = 2 4 / 3 ft,j^. 

In the frame of this model, in these two cases, the viscous force per unit length of meniscus does not depend on the 
meniscus radius r,„ [33] • 



B. Validity range of the models 



1. Velocity range 

At very small capillary numbers, or at a very small meniscus radius of curvature, the wetting film thickness predicted 
by eq. [8] becomes so small that the disjoining pressure can not be neglected anymore. Teletzke et at have shown 
in |54(, that the wetting film and meniscus profiles become independent of the capillary number for very small Ca, 
when they have the same shape as at equilibrium, as observed in [55j . We thus expect a viscous force that increases 
linearly with Ca, as is the case for non deformable interfaces in Stokes regime. 

At large velocity, inertia becomes important. Quere et al. [1% 5(| measured the thickness of the liquid film deposited 
on a fiber of radius rf pulled at constant velocity out of a liquid bath. They found a deviation from the theory for 
Ca > 10~ 2 . At larger velocity, they predict that the film thickness and the dynamical meniscus extension scale as 
hoo ~ r/Ca 2 / 3 /(l — We) and id ~ r/Ca 1 / 3 /(l — We), where We = pU 2 rf/j is the Weber number. The viscous 
force involves (in the stress free and rolling cases) the ratio id/h. As these two length scales have the same inertial 
correction, we thus expect that the viscous force departs from the theory at larger capillary numbers than the film 
thickness. Experimentally, the Ca 2 / 3 scaling has been observed for the viscous force up to We ~ 10 for films in 
cylindrical tubes [26| . 



2. Interfacial rheology 



The assumption of incompressible interface or of vanishing tangential stress are almost never perfectly verified and 
we estimate here the experimental conditions in which these conditions are achieved within the desired precision. 
The interfacial extension rate is given by 



dv(x, h(x)) 
dx 



X 



7 d 
Tj dx 



,d 3 h 
dx 3 



7 {ZpCa) 2 ' 3 W 



TTll 

f,OUt 



H 2 



(19) 



For the incompressible case, X — 1/2 and the extension is thus consistently 0. However, the extension reaches large 
values for the stress free case. The functions ^fH'/H 2 are plotted on Fig. [7] and have a maximum at 0.19 for the 
front meniscus and a minimum at -0.7 for the rear meniscus. For U = lem/s, r\ = 10~ 3 Pa.s, 7 = 30 10~ 3 N/m 
(i.e. Ca — 0.3 10~ 3 ) and r m = 2 10~ 4 m, the maximal extension or compression rate is thus of the order of 400 s _1 . 
The stress free model is thus applicable if no noticeable surface tension variations are observed for this extension 
rate. Considering an interface for which the interfacial stress variations are dominated by the interfacial dilatational 
viscosity rj s , and a desired precision of (a mt — j e q)/leq "C e, we obtain the validity criteria for the stress free model as 



-Ca 2j3 < 



(20) 



With the previous numerical values, the relative variation in interfacial stress is 1% for r\ s — 0.7 10~ 6 Pa-m-s. 

A similar approach can be used to determine the applicability range of the incompressible model. For the rolling or 
sliding cases, the problem is solved without specifying the tangential stress at the interface a mt . However, a mt can 
be determined a posteriori from the tangential stress continuity ([3|). Substituting the velocity field defined in eq. |4]in 
the relation [3] leads to 



da 1 



dx 



h d3h n 



nui 

X) - /i — - 

7 h 



7 (3/gCa) 1 / 3 
r H" 



(1-A) 



1-H 



3(3 H J 



As expected, a mt has a constant value for the stress free case (A = 1, /1 = 0). When A = 
variations of order Co 2 / 3 occur. Indeed, by integration over x between —00 and i, we get: 



Sa mt = 



= 7 (3/3Ca) 2/3 



(1 " ~ 3$ Gf/r{L) 



(21) 

1/2 (incompressible case), 
(22) 
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This incompressible model will be suitable for interfaces able to produce such interfacial stress variations without 
noticeable surface extension. In the rolling case (fj, = 0), the interfacial stress difference Sa mt equals the viscous force 
per unit length exerted by the wall. This is a direct consequence of the fact that the velocity field is a parabola 
centered in the middle of the film. For the same parameter values as above, we get Sa mt — 0.75mN/m. In the 
case of an interface dominated by elasticity, Sa mt — E5A/A, and for the interface extension to be less than 1%, the 
dilatational modulus of the interface, E, must be larger than 75 mN/m. More generally, in the rolling case, SA/A -C e 
if 

^Ca 2/3 <e, (23) 
E 

which is thus the validity criteria of the incompressible model, at the desired precision e. 
The sliding case involves the wetting film length and will be discussed in the next section. 

3. Meniscus size and length of the wetting film 

The dynamical meniscus is defined as the place where the viscous friction takes place (in the rolling or stress free 
cases). Figure [7] allows us to determine the extension t& along x, of this dynamical meniscus. Imposing the restriction 
that only 5% of the viscous force is localized on each side of the dynamical meniscus, we find that the dynamical 
meniscus, in dimensionless form, corresponds to the X — X p range [—32; 3] for the front meniscus and [—75; 3] for 
the r* one. So £ d = 35/ioo/(3/3Ca) 2 / 3 and l r d — 78/ioo/(3/3Ca) 2 / 3 . In physical units, for a typical capillary number of 
Ca = 0.3 10 -3 , this leads to an extension of the dynamical meniscus of the order of 0.2 r m , from the apparent contact 
point toward the wetting film. In a bubble geometry, as in Fig. [SJ the two menisci are thus independent if I > 0.2 r m 
at Ca = 0.3 10~ 3 . This condition becomes less restrictive for smaller Ca. 

On the meniscus side (X — X p < 0), we obtain an extension which is of the order of the meniscus itself for 
Ca = 0.3 10~ 3 . The fluid flow in the meniscus is thus probably not always n eglig ible. This may be why the viscous 
force measured for lamella is higher than the force found for isolated bubbles |25l l26l l28j . 

In the sliding case, not only the lower bound for the wetting film length should be considered, but also an upper 
bound. In this specific case, the interfacial stress variation (1221) increases rapidly with t and is dominated by the second 
term for long bubbles. The model derived from eq. Q]is not valid for large variations of the interfacial stress. Indeed, 
in such cases, the Laplace pressure derivative should involve the derivative of the curvature, but also the derivative 
of 7. However, an order of magnitude of the maximal length compatible with the sliding regime is estimated by 
assuming that 5a mt is of the order of j eq . Injecting this condition in eq. [pleads to: 

^raax ~ 2r m C(Z ^ (^4) 

For Ca = 0.3 10~ 3 this leads to £ ma x ~ 30r m . This theoretical limit has been observed experimentally in [57| . 

III. PRESSURE DROP FOR AN ISOLATED BUBBLE OR A LAMELLA 

A. Force balance 

1. Single bubble 

We now consider a single bubble in a capillary of constant section S and of perimeter V. We assume that the 
bubble is long enough to create a thin wetting film between the bubble and the wall (see Fig. The problem is not 
invariant in the z direction anymore. However the results of the previous section concern the flow in a small region 
between the meniscus and the wetting film of extension r m Ca 1 / 3 . These results thus remain valid as long as the tube 
dimensions are much lar ger than rVnCa 1 / 3 . Some corrections should nevertheless be considered if the tube has an 
angular cross-section [22l l37l. [H, HH . 

In Fig. [3J we define three planes perpendicular to the tube section: A up and Ad own are close to the rear and 
front meniscus, but outside the dynamical meniscus, and Ai is the plane in the middle of the bubble, situated at a 
distance I from each effective contact point between the meniscus and the wetting film. As already discussed for the 
film profiles in Fig. [5j if the bubble is long enough, the wetting film properties close to Ag are very close to their 
asymptotic values at +00, which are identical for the front and r* films. For menisci in a sliding regime, t represents 
the length of wetting film associated with each meniscus (see section III AI) . 



11 




Solution 



Tube 

FIG. 9: Sketch of a single bubble in a tube, with the two subsystems & and ST used for the force balance. The two subsystems 
are defined by the liquid bounded by the tube wall, the median plane At and, respectively, the plane Adomn and A up . 



The force balance made on the open system Qf (resp. f2 r ») constituted by the liquid between Ad 
and Ai, leads, if momentum fluxes are negligible, to the relation 

P f/r =P°-'^(7(i) + f f/r ) 



(resp. A up ) 



(25) 



with p° the pressure in the bubble and pf' r the average pressure, respectively in the section Ad OW n and A up . In the 
static case, f 3 ' l r = 0, and the relation predicts the Laplace under-pressure in the liquid phase. In the dynamical 
regime, the pressure in the liquid is smaller than the static value for the front meniscus and higher for the rear 
meniscus, as ff > and f r < (see eqs. pSlfTT])). Note that the pressure and velocity fields have been obtained in 
section [TT] by assuming a static pressure in the meniscus. Thus, obviously, the dynamical pressure drop can not be 
deduced from the pressure distribution plotted on Fig. [7J The dynamical pressure drop appears as a higher order 
correction that could in principle be re- injected as boundary condition for the film profile, in an iterative scheme [lij ]. 

The total pressure drop for an isolated bubble is Sp = p r — p* = (V /S)(ff — f r ) — (V /S)f '. For the case of a 
cylindrical tube of radius r t the previous relation simply becomes dp = 2/ /r t , so from eq. 1181 



Sp sf = 4.94 Ca 2/3 



27 



2"/ 



5p sl1 = 6.82Ca 2 / 3 + 0.94 — Ca 1 ' 

\ r m J r t 

Sp roL = 6.22 Ca 2/3 ^2 
n 



(26) 
(27) 
(28) 



The radius r t and r m are identical for the bubble, but differ for the lamella case, discussed below. 

This pressure drop, related to the presence of the meniscus, can be compared to the pressure drop induces by the 
Poiseuille flow in the tube, on both sides of the bubble Sp Pms = 8r]L t U/r 2 , with U the average liquid velocity and 
L t the tube length. For r t = 1mm and Ca = 10~ 3 , the meniscus contribution (in the stress free case) is of the same 
order than the pressure drop induced by a length L t = 12mm of tube. The local force induced by the meniscus thus 
strongly dominates the confined diphasic flows, even for relatively distant bubbles. 



As discussed in section Til Al the sign of the viscous force is not obvious and may vary for the rear meniscus. An 
energy balance is made below (for the stress free case), that provides a condition for this sign 
First, the mass balance on or fl r leads to 



Vh^U = SU' 

with U' the average velocity on the plane Ad OW n or A up . The energy balance is 



(29) 



V 



f/r SU' - p^Vh^U + Vf f/r U + V-yU - VjU > . 



(30) 



The last two terms are respectively the power provided by the interface and the free energy flux exiting the open 
system f2. Replacing eq. 1251 and eq. [551 into this equation leads to the condition, at dominant order 
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— >u 

FIG. 10: Sketch of a lamella. 



Uf f/r >U7 h^. (31) 

If U > 0, this leads to f* > r )h, x T' /S > 0. In contrast, for U < 0, we get f r < jhooV/S. The viscous force can thus, 
a priori, take positive or negative values. Using ft M = H"{— oo)r m (3Ca) 2 / 3 (eq. [7]) and f r — r y(3Ca) 2 ^ 3 F r (eq. [T3")) 
and taking V jS = l/r m which is the most restrictive condition, we get 

F r < H"(-oo) (32) 

which is verified by the numerical result shown on Fig. |SJ 



2. Lamella 



Considering the case of a single lamella, the situation only differs in a few details. The rear meniscus is downstream 
and the front meniscus upstream. A first consequence of this is that the wetting film downstream can be of arbitrary 
thickness, and the viscous force corresponding to the adapted rear profile must be used (see fig. |5]). A force balance 
carried out on the open system depicted on Fig. [10] leads to the relation Sp = Pq P — Po° wn = (2/ r *)Cf — / r + 7^ iy) ~ 
r y r {( r )) for a cylindrical tube. In contrast with the bubble case, the force balances made for the rear and front menisci 
involve the surface tension at two distinct points and the difference — r y r (£ r ) is thus a priori unknown for the 

incompressible cases. However, if the external boundary conditions impose the same surface tension at both sides of 
the lamella, we get the same results as for an isolated bubble, eqs. (|26H28I) . 

For lamella, another limitation must be considered, that the meniscus is not pushed by a pressure difference, but 
pulled by the lamella. The force exerted by the lamella, per unit length in the z direction and projected in the x 
direction, is 27 cos ip, with ip the apparent angle between the lamella and the wall, defined in Fig. [TUJ This force per 
unit length can not exceed 27, which thus determines the maximal velocity U m at which a lamella can be pulled 26]: 



f{U m )=2 1 . (33) 

In the stress free case, this leads to the simple condition Ca < 0.26, or, for 7 = 30 10 -3 N/m and 77 = 10~ 3 Pa-s, 
U m = 8m/s. Above this maximal velocity, the lamella breaks. For foam transport in porous media, as used for oil 
recovery or soil depollution, this dynamical criteria for the lamella stability may be an important issue. 



B. Bretherton's pressure derivation for a single bubble 



In the specific geometry of a single bubble in a cylindrical tube, the pressure can be obtained using a completely 
different approach, based on both the asymptotic interface profile shown in Fig. [4] and on the static Laplace relation. 
This has been derived by Bretherton [13| and in a more rigorous formalism by Park and Homsy flij ]. 

In the small capillary regime discussed in this paper, the bubble is almost at equilibrium far from the wall. Its 
shape is thus a spherical cap. As discussed in section llBl the extrapolation of the outside profile reaches its minimum 
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at the position h p = r m (3/3Ca) 2//3 H'j out H p (positive or negative). The actual radius of curvature of the meniscus is 
thus slightly smaller than r t at the front and slightly larger at the rear side. At dominant order in Ca, the dynamical 
radius of curvature is thus 

rU r ' = n (l + WCaf!*Hl out H^*) (34) 

The Laplace pressure becomes: 



2 7 (-oo) 2 



r f/r " n 



P°-- (7(-oo) + ^Caf^H'l^H^') . (35) 



This equation is the same as eq. [25] (with V jS = 2/r t ), but this equivalency requires some discussion. In the stress 
free case, 7 is everywhere equals to 7 e(? . We deduce from eq. [13] that eq. [25] and [35] are identical only if 

= H'j tOUt Hf' r . (36) 

This relation can easily be proved by integration by parts of F^ r , using the relation [5] 

For incompressible interfaces, the situation is more subtle, as the correction in the surface tension is of the same 
order as the correction in the curvature. Thus both need to be taken into account in dp. The precise value of j(— 00) 
is given by a 1 "^. Using, as previously, <j int 'f(t) — a %nt ' T (£.), Sp becomes 

Sp = ™ (So**'* 6a int > r ") ~ ^(3(3Ca) 2 / 3 H'l out (Hf - ff/) . (37) 

Inserting eq. [22] and [36] into this latest expression leads to the same result as given by eq. [27] and [28] thus insuring 
the consistency of both approaches. 

The pressure correction induced by the surface tension variation is not considered in [60j . thus explaining the 
different pressure value proposed for the rolling case in Fig. 6 of j60]. It should also be noted that using a direct 
analogy of this method to obtain the pressure drop for a lamella, as made in [l| , leads to a wrong relation between 
the friction force / and the meniscus radius. 



IV. BOUNDARY CONDITIONS FOR FOAM FLOWS 



A. Film orientation 



The friction between a liquid foam and a smooth solid wall is governed by the force exerted on each meniscus in 
contact with the solid. The model previously developed only considers a meniscus moving in a direction perpendicular 
to its axis. However, in a foam, the menisci in contact with the wall are oriented in all directions. Comparison of the 
pressure drops associated with lamellae of different orientations and moving in channel of rectangular section, suggests 
that only the projected length, in the direction perpendicular to the flow, contributes to the pressure drop [24], |6l|, |62j . 
However, this conjecture, that was natural for a foam having a motion of solid translation, is difficult to extrapolate 
in a general case. Indeed, during a complex deformation, only the normal velocity of the film is experimentally well 
defined. It is thus more consistent to assume, as in the viscous froth model [63|, [64[, that, locally, the force is oriented 
perpendicular to the meniscus and depends only on the normal velocity. We consider a portion of meniscus of length 
ds, that makes an angle 9 with the x direction and that moves with the apparent velocity Ux (see Fig. [TT1) . The 
viscous forces dF x projected in the x direction acting on this element are respectively, for the projected length model 
(pi) and for the non linear viscous force model (vfm), 

dF£ l = -fcos6ds, (38) 

and 

dF vfm = _ e ^£y\ cose yo/3 ds _ ( (^y 3 A (cos0 )4/3 dS; (39) 

with / the viscous force per unit length of meniscus as defined by ea llSI 

No clear experimental evidence allows us to discriminate between these two models yet. For the sake of simplicity, 
we will use eq. [35] in the following section. 
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FIG. 11: Network of menisci in contact with a wall, for a 2D or 3D dry foam. Each meniscus makes an angle 6 with the 
direction z, normal to the foam velocity. The contact area between a bubble and the wall is A c . 



B. Influence of the liquid fraction 

The effect of the liquid fraction is implicit in the sliding case, through the ratio £/r m (see eq. Additionally, the 
tangential stress between the foam and the wall depends on the total amount of contact lines between the meniscus 
and the wall, per unit surface, which strongly depends on the bubble volume V = 47rr^/3 and on the liquid fraction 
o 12*1. 



1 . 3D dry foams 

In the dry foam limit, the contact between a bubble and the wall is a polygon of area A c (see fig. [Til . The 
statistical relation between the bubble radius and its contact area has been investigated numerically in [65[, and, 
for a monodisperse foam, the average contact area is found to be A c = 2.5 r% . Assuming a hexagonal shape of arbitrary 
orientation for the contact area, we get the averaged projected length of meniscus per bubble, d pro j, as derived in 

/ 2A 

d pro] = 3(| cos0\) = 1.87r 6 (40) 

Here we neglect the area occupied by the meniscus itself on the solid wall. For slightly larger liquid fractions a 
correction has been proposed in 28] . 

If the interfaces are incompressible, the force / depends on r m and (£}, both of which must be determined as a 
function of 4> and TV In the limit of low liquid fraction, the radius of curvature of the meniscus is well approximated 
by r m = 1.73rt\/?> [13] ■ Computing the average distance between both sides of a regular hexagon (with 2 edges 
perpendicular to the velocity direction) we obtain (£} — 0.64rj, and 

(£)/r m = 0.37<T 1/2 . (41) 
Finally, for a monodisperse dry foam, the tangential stress exerted by the wall is 

a dry = %^/(Ca, (£)/r m ) = —f(Ca, O.370- 1 / 2 ) , (42) 

with the expression of the force per unit length of meniscus / given by eq. 1181 for each kind of interface, and the 
bubble radius. 



2. Wet foams 



In a wet foam, the bubbles are much less deformed and their shape remains closer to a sphere. The area per bubble 
on the wall is thus close to Af, — irr%. The contact area is a disc of area A c — irr^, as depicted on Fig. [12] The 
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FIG. 12: Wet foam in contact with a wall. Each bubble occupies a hexagonal domain of area At = i^r\. The contact area 
between the bubble and the wall is a disc of area A c — -kt\ represented in grey. 

fraction of the wall a 2 in contact with the bubbles has been measured by Princen [66[ : 

r 2 3 9 

a 2 = % = 1 — „ . (43) 

r 2 / 1 i \ 0.5 



The total amount of meniscus (with its rear and front part) per bubble, projected in the direction perpendicular to 
the motion, is thus given by d pro j — 2r c = 2rha and the average wetting film length by (£) = ^ar^. Finally, following 
the approach proposed by Denkov et al. in [291 ] , the radius of the meniscus r m is determined from its relation with 
the osmotic pressure of the foam H OS m [6?| : 

n osm = a 2 7 k , (44) 

with k the interface curvature, outside the contact area. In this wet regime the bubbles are close to spheres with flat 
patches, giving k ~ 2/r m . So we can express r m as r m — 2^ya 2 /H osm and 

W A™ = ^n osm . (45) 

The osmotic pressures measured by Princen and Kiss (68j are described by the phenomenological relation 

7 0.00819(1 - <j)) 2 



L osrn,i J r 



r b (0.0361 + O.96390) 2 



(46) 



Holder et al. provided more recent numerical and experimental data for II osm in [69j . approximated by the slightly 
different relation 

7 (0.26 -cj)) 2 

Ltosm,Ho — I O = . (47J 

These relations enable us to predict the tangential stress at the wall as a function of the bubble radius ri, and the 
liquid fraction cf> in a wet foam: 

a wet = %^/(Ca, (£)/r rn ) = —f{Ca, (£)/r m ) , (48) 
with a and (£) / r m expressed as function of 4> a- n d r b in eq. [43] and 1451 and / given by eq. 1181 

C. Comparison with some experimental results 

1. 3D Foams 



Foam rheological properties are usually investigated with a rheometer in cone/plane geometry. Using smooth plates, 
there is a velocity range in which the foam is not sheared in bulk and only slips on the bottom plate. The relation 
between the slip velocity and the stress can thus be deduced, as a function of the bubble size and liquid fraction. 
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For foaming solutions leading to stress free interfaces, the scaling a ~ Ca 2 / 3 is well verified [29j]. However, the 
theoretical law obtained from eq. 05] and eq. ITfflmust be corrected by a prefactor close to 4. Similarly high prefactors 
have also been observed for 2D foam 128] and lamellae I24L 12611. 



Experimental data have been obtained by Denkov et at for several foaming solutions of very high interfacial 
dilatational modulus [29| . The results of Section II C 21 show that the film profiles obtained with this kind of foaming 
solution are well fitted by the sliding model of interface. We thus compare the experimental stress with the prediction 
from eq. 1481 using the expression of / corresponding to the sliding interface (eq. [H 



- ^-7 

^wet.sli ^ 

n 



6.82 Ca 2/3 + 0.94— Ca 1 / 



(49) 



with Z = 2a /n. The experimental liquid fraction is 

\Pr 



0.1. 



The two models of osmotic pressure eq. 25] and eq. |4"T1 
lead respectively to (£)/r m )^ r — 0.318 and (£)/r m )"° = 0.5. Both predictions have been plotted on Fig. Q21 keeping 
Z as an adjustable parameter. A very good agreement is obtained with Z^ zt - Pr — 0.7 and Zf lt ' Ho = 0.5, whereas 
the theoretical prediction for the prefactor is Z theo = 0.3. Consistently with the observations at the single meniscus 
scale, the model slightly under-predicts the force. However, the variation with the capillary number is quantitatively 
reproduced by our model, without adjustable parameter on (£)/r m . 




Ca 

FIG. 13: Stress at the wall rescaled by j/ri as a function of the capillary number Ca for a 3D foam. Symbols: Experimental 
data obtained by Denkov et al. in [29]. Fig. 8c. The solid line is obtained from eq. [39]with Z = 0.7 and £/r m — 0.318, and the 
dot-dashed line with Z — 0.5 and £/r m = 0.5. The dashed lines are respectively power laws of exponent 1/3, 1/2 and 2/3. 

Another set of experimental data obtained by Marze et al., with a solution of Amilite GCK-12, tests the variation 
of the stress with the liquid fraction in the range tf> = 0.1 — 0.25 [3(|. The data, plotted as a function of the 
capillary number, are well fitted by power laws, with an exponent n that depends on the liquid fraction: n = 
(0.28; 0.44; 0.48; 0.66) for <j) = (0.1; 0.15; 0.20; 0.25). Qualitatively, the exponent varies continuously from 2/3 for a 
dynamic dominated by the meniscus (high liquid fraction) to 1/3 for a dynamic dominated by the wetting film 
(moderate liquid fraction), as expected for incompressible interfaces in the sliding regime. A quantitative comparison 
of these data with the prediction of eq. HH1 is shown in Fig. Q3] using the expression eq. H7] for the osmotic pressure, 
as it gives a much better prediction for high liquid fractions than eq. 1461 The agreement is very good for the highest 
liquid fractions but the 1/3 power law obtained for the smallest value of <fi is not reproduced. In that case, the best 
agreement is obtained for Z' %t = 2.0Z theo (using a single adjustable parameter for the four data series). 



2. 2D Foams 



The rheological properties of a foam strongly depends on its organization at the bubble scale and the need to 
observe this local structure is at the origin of the numerous experiments on 2D foams. The bubbles are organized in 
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FIG. 14: Stress at the wall, rescaled by 7/r(, i _p lotted as a function of the capillary number for a 3D foam of Amilite GCK-12. 
A, A. o, •: Experimental data from Fig. 8 in [30], with liquid fractions 0=0.10, 0.15, 0.20 and 0.25, respectively. The solid lines 
are predictions of eq. 1491 using the theoretical values of (£} /r m = 0.5; 0.25; 0.09 and 0.004, as obtained from eq. 1451 and eq. [47] 
and Z fit = 2 Z theo = 0.6; 0.44; 0.32 and 0.18. The dashed lines are power laws of exponent 1/3 and 2/3. 

a single monolayer that is confined cither between two glass plates, between the solution and a glass plate or between 
the solution and the atmosphere. In the first two cases, the viscous friction at the plate plays an important role in 
the dynamics and is often several order of magnitude higher than the internal viscous force. The stress at the plates 
has been measured by Raufaste et al. as a function of the capillary number, the meniscus size and the bubble area 
[28j |. All their experimental data are well described by the phenomenological law 



with A c the area per bubble. The velocity power law in 2/3 for the force per unit length of meniscus is not consistent 
with a dependency on the meniscus radius. However, fitting only the velocity dependency, they found / ~ Ca ' 58 
which is between the 1/3 and 2/3 exponents. A finite interfacial compressibility may thus be at the origin of this 
behavior. 



In this paper, we show that a part of the available experimental results on meniscus motion can be explained 
using very simple assumptions for the interfacial properties. We focus on both the viscous force acting on the moving 
meniscus and on the induced pressure drop. These are quantities of high practical interest for controlling or predicting 
bubbles or foam motion close to solid walls. In particular, we establish new predictions for these quantities in the limit 
of incompressible interfaces. However, wetting film thickness measurements clearly evidence intermediate regimes, 
between the stress- free interface and the incompressible one [42l l57l FzoI . [7l| , and the influence of these on the forces are 
still unclear. The surfactant concentration at the interface is governed by convective and diffusive processes, as well as 
adsorption and desorption rates. The interfacial stress is modified by the surfactant's intrinsic shear and dilatational 
viscosities and by the variation of the surface tension. More realistic models, taking into account some or all of these 
interfacial properties, are developed numerically [H l42l l49l. I52I l60l. l70l. [72T - t75l | . Nevertheless, the rheological properties 
of the interfaces are not very well characterized at the high extensional rates generated by the meniscus motion 
and it is often difficult to identify the dominating processes. Moreover, for incompressible or almost incompressible 
interfaces, the outer flow, far from the dynamical meniscus, may be strongly modified by the moving interface [50i . 




The relevant boundary conditions required to model the interface motion are thus difficult to identify. Different 
assumptions for the outer flow may lead, in the limit of high surface modulus, to either the rolling or the sliding 
limiting cases. Consequently, solving the hydrodynamical problem at a larger scale, and not only in the dynamical 
meniscus domain, is probably unavoidable. 




(50) 



V. 



CONCLUSION 
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